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Abstract 

Observations across many families of unconventional materials motivative the search 
for robust mechanisms producing linear in temperature d.c. resistivity. BKT quantum 
phase transitions are commonplace in holographic descriptions of finite density matter, 
separating critical and ordered phases. We show that at a holographic BKT critical 
point, if the unstable operator is coupled to the current via irrelevant operators, then a 
linear contribution to the resistivity is universally obtained. We also obtain broad power 
law tails in the optical conductivity, that shift spectral weight from the Drude peak as 
well as interband energy scales. We give a partial realization of this scenario using 
an Einstein-Maxwell-pseudoscalar bulk theory. The instability is a vectorial mode at 
nonzero wavevector, which is communicated to the homogeneous current via irrelevant 
coupling to an ionic lattice. 



1 Unconventional d.c. and optical conductivities 



It is an observed fact that several chemically diverse families of unconventional materials 
exhibit a low temperature resistivity that is linear in temperature, when tuned to what 
appear to be quantum critical points. Illustrative recent experimental data can be found 
for cuprates, pnictides, heavy fermions, ruthenates and organic superconductors - see [1] for 
an overview and references. While it is likely that there may be more than one explanation 
for this behavior, the universality of the temperature dependence of the resistivity seen 
in the phase diagrams of these materials strongly motivates the search for robust physical 
mechanisms that can reproduce the observations. 

In this paper we show, in the context of the holographic correspondence, that the struc- 
ture of superradiant instabilities of extremal black hole horizons leads universally to a linear 
in temperature resistivity when these systems are tuned to the quantum critical point me- 
diating the instability. This mechanism will be described in some generality. We go on to 
illustrate the process in a particular model that captures additional features of several of 
the unconventional materials of interest: the resistivity will be due to scattering off modes 
that are becoming unstable, and which are supported at nonzero momentum. One can take 
these modes to model the spin density and charge density wave instabilities appearing in 
the phase diagrams of [1]. 

That the resistivity is linear rather than, say, quadratic in temperature is only half of 
the mystery. Unlike conventional metals, many of the materials of interest fail to exhibit re- 
sistivity saturation as the temperature is increased. The resistivity increases unabated with 
temperature through the Mott-Ioffe-Regel limit. The materials are consequently known as 
bad metals [2]. This fact becomes particularly confusing when considered in conjunction 
with the dependence of the in-plane (frequency dependent) optical conductivity on temper- 
ature. In conventional metals, the optical conductivity exhibits a Drude peak that broadens 
as the temperature is increased. Resistivity saturation occurs when the spectral weight is 
smeared out over the whole of the bandwidth and the Drude peak effectively disappears 
[3]. In several classes of unconventional materials that do not exhibit saturation, the Drude 
peak is accompanied by an extended tail that falls off slowly at large frequencies up to the 
bandwidth scale [4]. At the would-be saturation temperature, the Drude peak is observed 
to melt into this broader feature [3, 4]. The d.c. resistivity, however, continues to increase 
linearly in temperature with the same slope, irrespectively of whether or not there is an 
associated Drude peak! The conundrum is served: Does the linear in temperature resistivity 
originate in Drude-like (i.e. momentum-relaxing) scattering or not? 
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The answer to this question is again likely to be material dependent. The picture of 
resistivity saturation outline above may not be correct. Nonetheless, the characterization 
of bad metals as metals that do not exhibit a zero frequency collective mode encourages 
the notion that the resistivity might be controlled by quantum critical physics, presumably 
responsible for the extended tail in the spectral density, rather than be sensitive to the 
mechanism of momentum relaxation. Such a picture is likely to have trouble with the fact 
that at lower temperatures a very sharp Drude peak is observed on top of the broader 
feature, see e.g. [5] for measurements in optimally doped YBCO. The mechanism of linear 
resistivity presented in this paper will be quantum critical in nature, and we will assume 
that the Drude peak has been swamped by the critical degrees of freedom. 

Figure 1 below illustrates the above discussion. 
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Figure 1: Theorist's schematic view of the optical conductivity in bad metals at lower 
temperatures (left) and higher temperatures (right). As the temperature is raised, spectral 
weight is shifted from the Drude peak into the broad tail and to interband energy scales. 
The linear in temperature d.c. resistivity does not notice the melting of the Drude peak. 



2 Extremal horizons and BKT transitions 

In quantum theories with a holographically dual description, the strongly interacting physics 
is described by classical gravitational dynamics in a dual spacetime with one extra spatial 
dimension. The extra dimension geometrically implements the renormalization group flow. 
In particular, the far interior of the spacetime describes the far IR, or lowest energy scales, 
of the dual quantum system [6]. In many holographic models placed at a finite charge 
density, (J*) 7^ 0, the far IR geometry is characterized by an emergent scaling symmetry 
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that indicates a power law specific heat at low temperatures, and hence the presence of 
gapless degrees of freedom [6]. Working in two spatial dimensions for concreteness, the 
general scale invariant form of the metric is [7] 

This spacetime admits the scaling symmetry {t,p} — > \{t,p},x — > \ l l z x. Therefore z is 
the dynamical critical exponent. For simplicity we will not discuss here the more general 
class of metrics describing hyperscaling violation [8, 9], although our considerations apply 
to those spacetimes also. Holographic IR scaling spacetimes typically arise over a finite 
range of parameter space and therefore describe critical phases. 

In a scaling geometry of the form (2.1), all the operators O in the low energy theory have 
an energy scaling dimension A. This dimension determines the spectral weight (imaginary 
part of the retarded Green's function) in the regime oj <C T <C n to be 

lim -ImGg (w,T) ~T 2A - 2 - 2 ^. (2.2) 

uj—>-0 OJ 

Here the chemical potential p, has been used to indicate the UV scale. Physics well below 
that scale is captured by the low energy spacetime (2.1) and is amenable to dimensional 
analysis. We also used the fact that the spectral weight is odd in frequency for bosonic 
operators. In general, the above spectral weight is evaluated at zero momentum, k = 0. A 
slight generalization is possible in the case of z = oo. In this case, space does not scale and 
so the momentum k is dimensionless and can be nonzero. The scaling dimension becomes 
momentum dependent: A(fc). 

According to the basic holographic dictionary [10, 11], each operator O is dual to a field 
4> in the IR spacetime. For simplicity, to start with, consider the case in which cf> is a scalar 
field with mass m that is not coupled to other fields at a linearized level. By solving the 
bulk wave equation in the geometry (2.1) and reading off the scaling behavior from the 
solution as p — > 0, one immediately finds 



2 + z 2 + z l 90 (2 + z) 1 . . 



The parameter v introduced here will appear repeatedly below. In the case of z = oo, the 
mass m 2 can be momentum dependent. In this expression we see that if the mass squared 
satisfies the generalized Breitenlohner-Freedman bound 

Lf> 2 >-^^, (2.4) 
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then the scaling dimension A is real. 

The quantity Lf R m 2 can be tuned by varying UV parameters. In particular, we can 
imagine tuning these parameters such that the mass squared drops below the bound (2.4), 
and the scaling dimension becomes complex. It is well understood by now that this triggers 
an interesting quantum phase transition in which the operator O condenses once the mass 
squared of <fi becomes too negative. We will return to the nature of the transition very 
shortly, but first we notice that precisely at the critical point, where the square root in (2.3) 
vanishes and v = 0, then from (2.2) and (2.3) we have 

lim-ImGg (w)~^. (2.5) 

This expression, which has been previously emphasized in [12, 13], will be at the heart of the 
linear resistivity above a quantum critical point that we will discuss in the following section. 
Previous interest in this expression is due to the fact that it resembles the spectral weight 
of the bosonic mode underlying the marginal Fermi liquid phenomenology of the cuprates 
[14]. It is not quite the same, however, because in general (2.5) only holds at k = 0. Even 
in the case of z = oo, because the dimension A is k dependent - this is why these systems 
were termed semi-locally critical in [15], rather than fully locally critical - the full spectral 
density will have a nontrivial k dependence, and (2.5) will only hold for one value of k at a 
time. The absence of a k dependence is crucial for the marginal Fermi liquid, as the mode 
is coupled to a Fermi surface, which has spectral weight at a nonzero k = kF that is set 
by UV dynamics. In contrast, our framework will operate entirely at the level of currents, 
which control the holographic charge dynamics at leading classical order in the bulk, and 
will not involve explicit discussion of Fermi surfaces. 

The association of complex IR scaling dimensions to instabilities was first made in the 
context of holographic superconductors [16, 17, 18]. In cases where the operator O carries 
a charge, the instability can be understood as a cousin of the superradiant instabilities of 
charged black holes, driven by pair production of quanta near the horizon and leading to 
a discharging of the black hole [6]. It was later realized that the quantum phase transition 
mediating this instability was of Berezinskii-Kosterlitz-Thouless (BKT) type. For instance, 
when the mass squared is just below the bound (2.4), the temperature below which the 
instability occurs scales like 

T c ~ fie-*/^ 2 . (2.6) 

Similar exponential hierarchies control quantities such as the condensate just below the 
critical mass squared. Such zero temperature BKT transitions were first discussed in [19]. 
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Unlike the conventional BKT transition, they can occur in any dimension and are not tied 
to an interpretation of vortex unbinding, but rather describe the merger of a UV and IR 
fixed point. These transitions were subsequently noted to be rather generic in holographic 
settings [20, 21] and to admit a 'semi-holographic' description [22, 12], in which the only 
role of holography is to provide a critical IR sector in which the transition occurs [23]. 
Our discussion in the following section will be essentially semi- holographic in nature, being 
independent of most details of the UV region of the bulk geometry. In §6.1 of the discussion 
section we will take a further step back from holography and comment on the validity of 
the result we have just described for general BKT transitions. 

There are two key points we would like the reader to take from the above. Firstly, 
that given a 'critical phase' with an IR scaling symmetry described by the metric (2.1) one 
can induce a quantum phase transition by tuning the dimension of an operator to become 
complex. Secondly, at the corresponding quantum critical point, the spectral density of 
this operator has the temperature dependence (2.5). This last statement holds at k = for 
finite z, and at some specific when z = oo. 

3 Mechanism of linear resistivity 

The d.c. conductivity is given by 

a = lim -ImG% JX (co,T) . (3.1) 

At a nonzero charge density (J*) and if momentum is conserved, this quantity is problematic 
because, in addition to the contribution (3.1), there is a delta function in the dissipative 
conductivity at oj = 0. In order to relax momentum over experimental timescales, the charge 
carriers must either be parametrically diluted or must interact with parametrically heavier 
degrees of freedom [24] . The result is the broadening of the delta function into a Drude peak. 
The contribution of the Drude peak to the d.c. conductivity is intimately connected to the 
mechanism by which momentum is relaxed. For instance, umklapp scattering in a Fermi 
liquid gives rise to the celebrated T 2 dependence of the resistivity. Instead, we would like the 
d.c. conductivity to be dominated by the universal quantum critical dynamics underlying 
the spectral weight (2.5). This can happen if the Drude peak contribution is swamped by 
the zero frequency limit of an extended tail that arises from scattering off critical modes. 
As we discussed in §1 above, this may be the case in bad metallic regimes. 

It is sometimes asserted that the presence of a Drude peak is synonymous with a quasi- 
particle description of transport. This is not quite correct. The essential requirement for 
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a Drude peak is a hierarchy of time scales, whereby the momentum relaxation timescale is 
much longer than any other timescale in the system. This is particularly clear in hydro- 
dynamic or memory matrix approaches, e.g. [25, 24]. In the presence of such a hierarchy, 
strongly correlated systems without a quasiparticle description will still exhibit a Drude 
peak. Conversely, the absence of a Drude peak simply requires that momentum is being 
dumped by the charge carriers at a rate comparable to all other interactions. In such cir- 
cumstances, the spectral weight from the delta function is transferred into the critical tail 
or indeed to interband energy scales. Strong interactions are presumably important here to 
maintain a metallic character and avoid localization [2]. In the concrete model we consider 
in the following section, we assume that such a process is occurring in our system, without 
disrupting the momentum-conserving interactions that we consider, so that in effect we can 
ignore the delta function contribution to the conductivity. More generally, we can imagine 
that momentum-relaxing processes are already part of the critical system that is undergoing 
the quantum critical BKT transition. 

If J x itself were the operator undergoing the BKT transition, then combining (2.5) and 
(3.1) would directly give a linear in temperature resistivity at the critical point. Such a 
phase transition would correspond to the spontaneous generation of a uniform current and 
presumably requires spontaneous symmetry breaking of the global U(l) symmetry. We are 
not aware of holographic, or other, models where this occurs. However, it is easy for the IR 
critical behavior (2.5) to get communicated to the current operator via operators that are 
irrelevant from the IR quantum critical point of view. 

The IR scaling geometry (2.1) is generically deformed by irrelevant operators that drive 
a renormalization group flow up towards the finite density UV fixed point or cutoff. In the 
IR scaling regime, before these irrelevant operators kick in, we can imagine diagonalizing 
the equations of motion for a generic perturbation of the background to obtain decoupled 
gauge invariant fields Near the boundary (p — > 0) of the IR geometry, these satisfy 

d> J (p) ~ c'(> +2 /*- A ' + p A ' g?(u>,T)) . (3.2) 

Here Gf(uj, T) is the IR Green's function, the c 1 are constants, and Aj is the IR dimension 
of the operator. We have allowed a temperature T <C n that only affects the IR spacetime. 
One of these dimensions Aj is assumed to undergo a BKT transition of the form described 
in the previous section. 

Away from the IR geometry, the irrelevant operators will typically couple these pertur- 
bations. However, in the regime of interest oj, T <C /x, this coupling does not introduce any 
additional non-analytic temperature or frequency dependence. In appendix A we show that 
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a generalization of the usual matching procedure of e.g. [26] implies that to leading order 
at low frequencies and temperatures 1 

TmG%j*(u,T) = ^d*Tmg?(u,T). (3.3) 
/ 

Here the real coefficients d 1 will generically be nonzero if the irrelevant operators mix the IR 
mode with the gauge field A x . This is by no means automatic; the model of the following 
sections will achieve the required mixing by combining several interesting ingredients. In 
more familiar condensed matter language, we can think of the direct coupling in (3.3) 
between the current and a fluctuating order parameter as a cousin of the Aslamazov-Larkin 
process [27], in this case mediated by irrelevant operators. The key fact is that we couple the 
fluctuating field directly to the current, not going via e.g. a fermionic self-energy. It is now 
immediate that if one of these IR operators undergoes a BKT transition, then the spectral 
weight (2.5), plugged into the matching formula (3.3) leads to a linear in temperature 
resistivity at the critical point 

r = -~T. (3.4) 
a 

Note that, according to (2.2) and (2.3), the contribution from a general IR operator Oi to 
the conductivity is 

a = T- 1+2u < . (3.5) 

Recall that v was defined in (2.3). Assuming we are in a stable phase, then we see that the 

contribution of all the other operators with vj > are subleading in the sum (3.3) compared 

to the critical operator with v = 0. The above discussion is depicted in figure 2 below. 

While (3.5) allows the critical operator to dominate at the critical point, away from the 

critical point, on the disordered side, it implies that the exponent of the resistivity will be 

1 — 1v < 1. Here v is the exponent of the critical operator away from the critical point. This 

is in contrast to experiments in all the unconventional materials of interest, which show that 

when detuned from criticality, the exponent of the resistivity increases towards the Fermi 

liquid T 2 behavior [1]. To capture this behavior we would need to increase the scope of our 

model. The simplest way to do this would be to combine the critical tail contribution we 

are discussing with a more conventional Femi liquid Drude peak contribution. It may be 

that by coupling the critical operator to the Fermi liquid, along the marginal Fermi liquid 

lr To be precise, as explained in the appendix, equation (3.3) holds when the fi's are real, that is, on the 
stable side of the quantum critical point. When some of the vi's are imaginary, the general expression has 
a complicated logarithmic dependence on T and u. 



7 



r~T 




T c ~ |j exp(-Tr/V-v 2 ) 



T 



<0>* 





Quantum 
Critical Point 



Critical phase 



Figure 2: Schematic phase diagram. The BKT quantum phase transition occurs at the 
boundary of a quantum critical phase when the scaling dimension of an operator becomes 
complex, signaling a condensation instability. If the unstable operator is coupled to the cur- 
rent via irrelevant operators, then above the critical point the quantum critical contribution 
to the resistivity is linear in temperature. Close to the critical point, the mode becoming 
unstable has a strong effect on the d.c. and optical conductivities. 

lines of [12, 13] or otherwise, one can remove the T 2 contribution to the resistivity in the 
critical region of the phase diagram. For instance, very schematically, a form like 



would start to look closer to the data. 

Complications with the Drude peak are removed if one looks at the optical conductivity. 
The scaling arguments above now imply that the dissipative conductivity scales like 



a ~ vT~ 2 + T~ 1+2 



(3.6) 



-l+2i/ 



(3.7) 



This gives the marginal Fermi liquid form 



1 

1 



(3.8) 
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at the critical point. This may be compatible with observations in e.g. YBCO [14] and LSCO 
[28]. We discuss the optical conductivity in more detail below. It would be interesting to 
measure in detail the doping dependence of the power law tail of the optical conductivity 
in these materials, to assess the plausibility of a dependence like (3.7) close to the critical 
point. As we noted above, the expressions (3.7) and (3.5) are strictly applicable only on 
the v > side of the quantum phase transition. It is a feature of these models that the 
d.c. and optical conductivities have asymmetric behavior on opposite sides of the quantum 
critical point. 

There may seem to be a tension in the fact that J x does not acquire a vacuum expec- 
tation value and yet its correlator couples at the lowest frequencies to the unstable mode 
according to (3.3). We will verify explicitly in our concrete model below that these two 
statements are compatible. 

The mechanism of linear in temperature resistivity we have described is universal in the 
sense that it only depends on the onset of a holographic BKT quantum phase transition at 
the boundary of a quantum critical phase, combined with the presence of irrelevant operators 
that couple the IR critical operator to the current. The mechanism is independent of the 
details of the theory undergoing the transition and also of the UV completion of the critical 
theory. Beyond the tuning to the quantum critical point, no additional specification of 
dimensions of operators or dynamical critical exponents is necessary. In the remainder of 
this paper we describe a specific holographic realization of the scenario that we have just 
outlined. 

4 Lattices and finite wavevector instabilities 

A concrete holographic model realizing the mechanism outlined above can be achieved by 
combining three interesting ingredients that have been the focus of recent discussion: 

• Vectorial instabilities occurring at finite wavevector A;* > [29, 30, 31]. 

• A (semi-) locally quantum critical sector in which all momenta are critical [32, 33, 15]. 

• A lattice that is irrelevant in the IR scaling regime [24, 34, 35]. 

The essential idea is that the irrelevant coupling to the lattice will communicate the finite 
wavevector instability to the k = electrical current. The instability at /c* can have critical 
scaling because of the local quantum criticality. 
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The most important effect of a lattice is to broaden the delta function in the conductivity 
[24] into a Drude peak. Impurities will also achieve this effect [36]. Nonetheless, many of 
the interesting materials are believed to be very clean, and indeed, the Fermi liquid T 2 
resistivity observed over large ranges of temperatures in these materials, away from the 
critical regions, presumably originates from umklapp scattering off the lattice rather than 
impurities plus interactions [37]. As we have stressed repeatedly, however, in this work we 
are not interested in the Drude peak contribution to the conductivity. The role of the lattice 
for us will be to mix modes with different wavevectors. 

We will consider the following 3+1 dimensional bulk model [30] 

C =\R * 1 - \ * dip A dp - V (if) * 1 - \ t (if) F A *F - \ (cp) F A F , (4.1) 

where F = dA. This Lagrangian describes Einstein-Maxwell theory coupled to a pseu- 
doscalar field with both dilatonic and axionic couplings to the field strength. The corre- 
sponding equations of motion are given by 

iV = dn<pd*<P + g„uV-T Q<7^ F\ p F Xp - F^F v ^j , 
d(T*F + $F) = 0, 

d * dip + V' * 1 + -t' F A *F + F A F = . (4.2) 

We will assume that the three functions V, r and # have the following expansions 

K = -6 + imV + ... ) T = \-^ip 2 + #= + (4.3) 

It is sufficient to know the action for the pseudoscalar to quadratic order, as it will only 
appear as a perturbation of the background. In (4.3) we have furthermore set the asymptotic 
AdS^ radius to I? = 1/2. In the Lagrangian (4.1) we already set both the Newton and 
Maxwell constants to unity. All of these quantities can be scaled out of the equations of 
motion in this model. 

The key feature of the theory (4.1), to be described immediately below, is that it can 
develop vectorial instabilities, which involve the electric current. Closely related to this fact 
is that the instabilities of interest occur at nonzero wavevector A;*, avoiding the spontaneous 
generation of a homogeneous current, as we would expect. By subsequently introducing a 
lattice, that becomes irrelevant in the IR and does not disrupt the critical phase, we will 
communicate the effects of this instability to the spectral weight of the current operator, 
leading to strong temperature and frequency dependence in the conductivity. 
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4.1 IR spectrum of excitations and instability 



The equations of motion (4.2) admit the following AdS2 x R 2 black hole solution 



1 / dr 2 \ 
ds\ = -f dt 2 + — [— + dx 2 + dy 2 ) , A = (r-r+)dt, <p = 0. 




(4.4) 



Here the emblackening factor 



f = r 




(4.5) 



The horizon of the black hole (4.4) is at r = r + and the temperature is T = y/3r + /ir. The 
boundary of AdS% is at r — > oo. This will serve as our near horizon scaling geometry, as 
per our discussion in §2 above, which we have further placed at a finite temperature. That 
the effects of the temperature are restricted to the IR scaling geometry implies that T <C fi, 
the UV scale. The scaling regime described by AdSi x R 2 has dynamical critical exponent 
z = oo and an associated ground state entropy density. This ground state entropy has 
tainted the reputation of AdS% x M 2 as a ubiquitous tool in the applied holography effort. 
(Semi) local quantum criticality with z = oo is in fact compatible with a vanishing ground 
state entropy density if hyperscaling is violated [38]. Holographic scaling geometries with 
z = oo are the only currently known holographic duals (away from the probe limit) that 
share at leading bulk classical order the basic property of Fermi liquids of having spectral 
weight at low energies but finite momentum. It seems conceivable that AdS2 x M 2 may have 
the last laugh. 

As usual, perturbations about the background can be decomposed into transverse and 
longitudinal sectors that decouple from each other. The finite wavevector instability occurs 
in the transverse channel. We can write the transverse perturbations around the exact 
solution (4.4) as 



Sgty = h ty (t, r) sin (kx) , 
5 A y = a (t, r) sin (kx) , 



Sg X y = h xy (t, r) cos (kx) , 
6 if = w (t, r) cos (kx) . 



(4.6) 



The equations of motion (4.2) then yield the system of linear coupled equations 



-kdth X y + k 2 h ty - f (2 d r a + d 2 h ty ) = , 
2 d t a + I2kf d r h xy + d t d r h ty = , 
-r 1 d 2 h xy + 12 d r (fd r h xy ) + kf- 1 d t hty = , 
- f' 1 d 2 a + 12 d r (fd r a) - 12k 2 a + c x k w + 12 d r h ty = 0, 
r 1 d 2 w + 12 d r (fd r w) - (12k 2 + m 2 s + n) w + 12c±k a = . 



(4.7) 
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Introducing the new field <p xy through 



dt<t> X y — f d r h xy , 

the system of equations (4.6) leads to the linear system of equations 

(□ 2 -M 2 ) v = 0, 



(4.8) 



(4.9) 



with v = ((frxy, a, w) and the mass matrix 



/ 12 k 2 



M 



2k 



\ 



(4.10) 





144 /c 24 + 12& 2 - Cl k 
\ -Ylcxk 12k 2 + m 2 s +n J 

Note that the Maxwell field fluctuation decouples from the remainder at zero wavenumber, 
k = 0. The Laplacian appearing in equation (4.9) is with respect to the two dimensional 
metric 



ds 2 2 = -fdt 2 + 



dr 2 
12/ 



(4.11) 



The linearly coupled system of equations (4.9) can be diagonalized to yield three inde- 
pendent modes 



(□ 2 - 3) m = o 



dhi + y d r gi 



Mi )9i = , 



(4.12) 



12 f 2 * ' 12/' 

with //? being the three eigenvalues of the mass matrix (4.10). If we go to frequency space by 
writing gi(t,r) = e~ twt Ui{r), the solution with infalling boundary conditions at the horizon 
[39, 40] is then 



Ui{r) 



r+J I>,)r(l + ^)/(r) •^r +r -^ 2 F 1 ^,^± 1 



_+ 

'ii ^2 



I - V) . — r | — (z/j O -Vi) , 

(4-13) 

where a\ = \— 2^7=^ v% and = \\Jl + By expanding the above solution near the 

boundary of the AdS 2 , i.e. as r — > 00, we can read off the AdS 2 retarded Green's functions 
in the usual way (e.g. [26]) to obtain 2 



2/A 



r (1 - Vi ) r(±-i^ T + Vi 
r (1 + v^ r (i - 1^ - vi 



(4.14) 



The case Vi = 0, which will be of particular interest below, should properly be treated independently, 
with the solution to the wave equation expressed in terms of a Legendre function. Correct result are obtained 
by continuation of the general results to Vi — > 0. In particular, the spectra density lmG R ~ co/T at Ui = 0. 
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These have the expected form that is determined by the SL(2, R) symmetry of AdS2- Again, 
see for instance [26]. These are the IR Green's functions that will appear in the matching 
formula (3.3). The perturbation includes the transverse current mode 5A y (x), at nonzero 
momentum k. Below we will couple this mode to the homogeneous current using a lattice. 
For small uj we have 



/ tt \^ r(i-^)r(i + ^) 



\2^fl j V {I + Ui)T {\ - Vi ) 

+ i ^-.r(j + *)r(i- 

2 Uvs) r(i-^)r(i + ^ ( ) v ' 

If Ui is real, then taking the imaginary part we recover the scaling with temperature that 
we anticipated in (2.2) above, using the expression (2.3) for A in terms of v. In general, 
the eigenvalues of the matrix (4.10) are slightly complicated functions of the wavenumber 
k, but are easily found numerically. 

To illustrate the features of the spectrum of our theory, consider the special case where 
nig + n = 0. (We will in fact consider a different case for the numerics below.) In that case 
we have 

Ul =1^/1 + 4*2, 



^ =1 W 5 + 4 fc2 - ^ ^12 + (24 + cf) , 



U3 = 2 V 5 + 4 k2 + v 12 + ( 24 + c i) k2 • (4 - 16) 



From the expressions above we can see that v\ has two minima at 

1 Cl ^48 + cj 
'2^12 ^J2A + cj ' 



(4,7) 



at which 



^2min = 7A/ 12 " -T-^rr^- ( 4 - 18 ) 



192 

3 24 + c 2 

For ci > 2^/6 we see that there is a range of A; for which z^2 is imaginary. These are the finite 
wavenumber instabilities that we shall use. For our purposes, it does not matter whether 
or not the range of unstable momenta extends down to k = 0. What is important is that 
there are finite wavenumber instabilities involving the Maxwell field SA y . 

Even when the system is stable, when < c\ < 2\/6, we see that, as opposed to the 
pure Einstein-Maxwell case, the dominant mode will have ^2 min < 1/2. This renders the 
finite wavenumber operators relevant in the IR. A lattice of such operators would lead to 
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a strong backreaction on the IR geometry [24]. In the following subsection we will instead 
introduce a lattice for the charge density J*, that will be seen to be irrelevant in the IR. 

Within the IR geometry, neither c\ nor the unstable range of k is tunable. The wavenum- 
ber k is dimensionless in the IR because of the local criticality. However, once the AdS% x M 2 
geometry is realized as the IR of a full asymptotically AdS^ spacetime, then we shall see 
that k must appear in the UV-dimensionless combination fe/\/j*, with J* being the total 
charge density in the system. By varying the charge density, which can be thought of as a 
proxy for varying the doping of an experimental system, we effectively vary the momentum 
that appears in the scaling IR Green's functions. In the following subsection we shall intro- 
duce a lattice wavenumber k^. By keeping the lattice fixed and varying the charge density 
J*, we will induce a BKT quantum phase transition in modes that couple via scattering off 
the lattice to the homogeneous current. 

The scenario outlined in the previous paragraph, and to be fleshed out below, does not 
quite give a realization of the phase diagram of figure 2: the system is either stable or 
unstable over some range of momentum and this fact cannot be tuned. What we will vary 
by varying k/VJ* is whether the instability feeds through to the homogeneous current, our 
observable of interest. To truly realize the desired phase diagram of figure 2 we need to 
enlarge the model to allow e.g. the parameter in the action c\ or the radius of the IR AdSi 
metric to be tunable from the UV. We do not see an obstacle to doing this, but will continue 
with our slightly simpler model in this paper. 

4.2 Lattice 

The previous subsection studied perturbations of a finite temperature scaling solution to the 
theory (4.1). Following [30] we exhibited an instability of this background, for c\ sufficiently 
large, over a range of momenta, kj± < k < ks- In this range the exponent vi (k) in (4.16) is 
imaginary. At the boundary of this range of momenta viO&aib) = 0) leading to the universal 
spectral function (2.5). In this section we will couple this critical mode to the k = current 
by introducing a lattice that is irrelevant in the IR. We will furthermore explain how k can 
be tuned to the critical values by varying the charge doping of the system. 

The simplest and perhaps most natural lattice to consider is an 'ionic lattice' in which 
the operator that couples to the lattice is the charge density J*. That is, we introduce a 
spatial modulation to the chemical potential 

Af ] =fi + V(x,y). (4.19) 
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Holographically this is implemented by the usual UV boundary condition lim^oo A t = A\ . 
In a purely electric background, we can consistently set the pseudoscalar ip = in the theory 
(4.1). Therefore the background generated by the lattice is found by solving the Einstein- 
Maxwell equations (sans pseudoscalar) subject to the UV boundary condition. We shall 
do this numerically in a perturbative approximation in the following section. However, the 
essential physical features can be determined without an explicit background solution, as 
we proceed to show. 

The charge density at any nonzero wavevector is an irrelevant perturbation about the 
AdS2 x M 2 solution (4.4) of Einstein-Maxwell theory [41, 24]. Therefore, however strong 
the UV potential in (4.19), the effects of the lattice will always be small in the far IR. The 
IR Green's functions computed in the previous section are thus unchanged. The role of the 
lattice is only to mix the unstable finite wavenumber modes with the homogeneous current 
mode, according to our discussion in §3. To isolate this effect from Drude momentum- 
relaxation physics, we take the current to propagate in the y direction, but the lattice to 
oscillate only in the x direction, with wavenumber k^ 



According to (3.1), to obtain the conductivity we must calculate ImGjy Jv (uj , T) at 
k = 0. We are now taking the current to run in the y rather than x direction. This Green's 
function is holographically related to perturbations SA y in the bulk. Even without a lattice, 
in a nonzero charge density spacetime, 5A y couples to Sgty With the lattice, inspection 
of the equations reveals that in the full spacetime the four modes {5A y , 5gt y , Sg xy ,S(p} are 
coupled at wavenumbers that are integer multiplies of k^. These can be collected into three 
gauge-invariant modes for each wavenumber 



satisfying second order equations of motion. In the IR, the different wavenumber modes 
decouple and obey the equations that we have solved in the previous subsetion. The upshot 
is that the general matching formula (3.3) becomes 



Here the IR Green's functions Qf- are those obtained in (4.14) above. Their wavenumber 
dependence is through Vi{k n ) in (4.16). The mixing is illustrated in figure 3 below. 




(4.20) 



k n — n kj^ 



(4.21) 




(4.22) 
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Figure 3: The IR is tuned to the boundary of an instability condensing finite wavenumber 
vectorial modes. The lattice is imposed in the UV but irrelevant in the IR. Away from the 
locally critical IR region, the lattice mixes modes of different wavenumber and couples the 
unstable mode to the homogeneous electric current. 

To tune the system across the instability identified in the previous subsection, insofar 
as the conductivity is concerned (see the discussion at the end of §4.1), we now need to 
vary k^ so that one of the k n becomes equal to the wavenumbers k^ or ks bounding the 
unstable region. Physically it is presumably not feasible to tune ki itself, as the material 
is fixed. However, in the formulae for Vi(k n ) in (4.16), a dimensionful quantity has been 
suppressed and lengths have been rescaled relative to natural lengths of the UV geometry 
[6]. This is possible because of the emergent z = oo scaling, under which lengths are 
dimensionless. The invariant dimensionless quantity to consider, from a UV perspective, is 
the ratio of momentum to the square root of the charge density: k/V~J^ [24]. The charge 
density J* corresponds to the total electric charge available for conduction. This can be 
tuned experimentally by doping the system. Doping is indeed the mechanism of tuning in 
the cuprate and pnictide unconventional superconductors. Thus, with kp fixed, doping the 
system will allow the dimensionless IR wavenumbers k n to be tuned through their critical 
values, communicating the BKT transition of interest into a linear in temperature resistivity 
(3.4). 

The optical conductivities in (3.7) and (3.8) now show that a lattice combined with 
local criticality and a mode on the verge of a finite wavenumber instability leads to a 
significant low energy spectral weight, ranging to a ~ a; -1 at the critical point itself. A 
strong optical conductivity going like a ~ a;~ 2 / 3 was recently obtained in an impressive 
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numerical holographic lattice computation in [34] . The mechanism underlying these effects 
are quite different. The lattice scatterings captured in [34] relaxed the current, resolving 
the Drude peak, while ours merely mix modes of differing transverse momenta. The tail in 
the optical conductivity observed in [34] is indeed a tail of the Drude peak itself, and would 
presumably not be there if momentum were conserved. Our tail is an additional feature 
due to scattering off IR critical modes. Remarkably the scaling a ~ cj -2 / 3 agrees with a 
reported measurement in BSCYCO [42]. As we noted above, a a ~ W 1 scaling has been 
suggested for e.g. YBCO [14] and LSCO [28], although a sophisticated scaling analysis of 
the data comparable to that in BSCYCO [42] remains to be performed. To substantiate 
these similarities, a connection to the underlying dynamics of the material is required. It is 
certainly of interest to combine these two effects of the lattice in a single model. 

5 Numerics: lattice deformed RN black hole 

The explicit realization of the mechanism for linear in temperature resistivity, and strong 
optical conductivities, of the previous section depended on various genericity assumptions 
about modes coupling to other modes. In this section we perform a numerical analysis with 
a lattice that will confirm all the assumptions made. To simplify the analysis, we work 
perturbatively (to second order) in the strength of the lattice. This allows us to use a large 
system of coupled ODEs rather than PDEs. The perturbative treatment is consistent and 
already contains all of the physics we wish to illustrate. 

5.1 The background lattice perturbation 

In the absence of a lattice, the full asymptotically AdS± background of the theory (4.1) 
dually describing a finite charge density, zero temperature state of matter is the extremal 
Reissner-Nordstrom black brane [6] 



We wish to deform by introducing a chemical potential or "ionic" lattice. We achieve this 
by setting the boundary value of the Maxwell field to be (cf. (4.20) above) 



ds\ 



f dt 2 + f~ l dr 2 + r 2 (dx 2 + dy 2 ) , 




(5.1) 



lim A t = 1 + A cos (Ulx) . 



(5.2) 
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In this subsection we set up the equations describing the effect of this deformation on the 
static background (5.1) to order A 2 . 

Choosing to work in a radial gauge where g rfl = 0, we can see that the components 
of the metric that will be affected by the lattice deformation are {gtt,9rr,9xx,9yy}- The 
nontrivial functional dependence of the metric components will be on the radial coordinate 
r and the spatial coordinate x. There is a residual gauge symmetry generated by the vector 
£ = £ r d r + £ x d x with V r £,x + Va;^ r =0. These perturbations of the background remain 
within Einstein-Maxwell theory (that is, the pseudoscalar is not sourced). As shown in [24], 
the u = and finite k perturbations of the near horizon AdS2 x M 2 geometry in Einstein- 
Maxwell theory are all irrelevant. This fact strongly suggests that the AdS2 x M 2 IR limit 
is not modified by the inclusion of a lattice. The magnitude of the effects of the lattice are 
tunable at the boundary and become small as r — > r + in the interior. We therefore conclude 
that as far as the background is concerned, for small parameter A we can consistently treat 
the backreaction of the metric and gauge field perturbatively at all radii. 

Expanding our field components in A we can write 

A t = 1 - r -± + A cos (k L x) + A 2 (a^ + cos (2ife L x)) + O (A 3 ) , 
9tt = -f [l + Aq£' 1} cos (k L x) + A 2 (q[ 2 ' 0) + Q (2 ' 2) cos (2k L x))] + O (A 3 ) , 



g rr = r 1 [l + A q£' 1} cos {k L x) + A 2 (qJ 2 ' 0) + q[ 2 ' 2) cos (2fc L x))] + O (A 3 ) , 

gu = r 2 [l + \Q$ 1) cos (k L x) + A 2 (Qf Q) + cos (2fe L x))] + O (A 3 ) , (5.3) 

where / is as in (5.1) and the a and Q functions depend on r only. We have used the 
residual gauge freedom mentioned above to fix the functions appearing in g rr in terms of 
the functions appearing in g tt . From the expansion above we have, computing the surface 
gravity, that the temperature is still given by 3 AirT = f'\ r=r . Plugging the ansatz for the 
perturbations (5.3) in the equations of motion we obtain 

• At order A: second order differential equations for the perturbations 
and an algebraic equation for QH' 1 ^ ■ 

• At order A 2 : second order equations for j<2^'°\ Q^' 2 \ a^ 2 '°\ a^ 2 ' 2 ^|, a first order 
equation for Q# and an algebraic equation for Q[ 2 ' 2 ^ ■ 

We will solve these equations numerically below. At this early point we note that for our 
purposes we will not need to solve for the functions Q^ 2 \ where a = t,x,y, and a^ 2,2 ) 



3 We will later impose that the metric perturbations fall off at infinity, so that the normalization of time 
remains the same as for the A = case. 
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but it should be noted that they are non-trivial by consistency of the equations of motion. 
In order to fully define the boundary value problem we now need to impose boundary 
conditions which will give us a regular near horizon geometry and with the right sources at 
the boundary of AdS^. 

The horizon is a singular point of the equations of motion and regularity of the solution 
imposes the following near horizon expansions at r = r + 

= QaaW + QaS) (r - r + ) + ■ ■ ■ , a = t,x,y, 
a™ = a f Q f (r - r+) + a^f (r - r+f + ■■■ . (5.4) 

Plugging the near horizon expansion (5.4) into the equations of motion, we discover that 
the constants , Q^gj an d Q^q) are cons t an ts of integration as well as a combination 
of Q^'Jj which we take to be the sum = + Q^y 

On the other hand, at infinity we have the following schematic expansions 

n (i,j) _ AhJ) , a (*>•?') r -i j 

n(i.i) - (7 {1 ' 1} r + G (1,1) + • • • + r~ 3 + • • • 

Virx - ^(-1) r + ^11(0) + + Lr xx(3) r + 

n(l.l) - r + <7 (1;1) + • • • + 5 r~ 3 + • • • 

^yy ~ U xx(-1) r + U yy(0) + + ^(3) r + 

<3?"=0§g + -. (5.5) 

In order to source the gauge field corresponding to the ionic lattice deformation (5.2) we 
need to set = 5 tl 5^ 1 A, which uniquely fixes the sources of the gauge field. The ionic 

lattice is the only deformation of the Reissner-Nordstrom black hole we would like to turn 
on and therefore we should also set G^f_ 1} = G-^oj = G fi(o) = G tt{o) = °- 

To conclude this section we note that we will have a total of six second order differential 
equations 4 in ^Q^' 1 \Q^ ,0 \a^ 1 ' 1 \a^ 2 ' ^ and one first order differential equation in 
We have eight constants of integration from the horizon {Q^gj, s^ 1 ' 1 ^, Q^q), a (o)^> a (0)°^} 
and five constants of integration from the asymptotic data j G^'j^ , G^ j , ^ji) 1 ^) } ■ 
This combination of differential order and number of constants of integration gives a unique 
solution for any given temperature. We will describe the numerical shooting method we are 
using in §5.3. 

4 As mentioned before, we don't need to solve for the functions Qia' 2 ' and a' 2 ' 2 '. 
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5.2 The A 2 corrected current-current correlator 



Given the background whose construction we have just described, we now wish to compute 
the retarded Green's function Gj y j y (uj,T) at k = in the deformed background. As will 
immediately become clear, for a source for J y at k = to induce a response also in J y at 
k = via scattering off a finite wavenumber mode, it is necessary to go to second order in 
the lattice strength, which is why we have also considered the background to order A 2 . To 
compute the Green's function, we consider the following time dependent perturbation 



5 A y = e 



Sip = e 
Sgty = e 
5g xy 



-iu>t 



-lUlt 



a(°'°) + A a*, 1 ' 1 ) cos (k L x) + A 2 (a™ + a y 2 ^ cos (2k L x)) 

sin (2kLx) + O (A 3 ) , 



+ O (A 3 ) , 



-iuit r 2 



dlT + cos (k L x) + A 2 (g^ + g^ cos (2k L y))] + O (A 3 ) , 

A ggjjV sin (k L x) + A 2 5 g 2) sin (2k L x)] + O (A 3 ) , (5.6) 



with all the functions 3>( l >.?') depending on the radius r only. The perturbation described by 
(5.6) is consistent and we chose again to work in radial gauge. We will fix the constants 
of integration in this gauge and then perform a large gauge transformation that will bring 
the perturbation into a form where it will be manifest that we are only sourcing the k = 
component of the current J y . 

Plugging the perturbation (5.6) into the equations of motion we obtain 

• Second order equations for jaj^, V^\9xy^} • 

• First order equations for {<?ty^} • 

Regular, infalling boundary conditions at the horizon gives the near horizon expansion 





= r 


4ttT 








op 


= r 


4irT 


yxy 


= r 


4ttT 


(hj 

Sly 


= r 


iuj 
4ttT 



(i.j) , (i.j) 



a y(0) + a y(l) 



(r - r+) + • • 
(r - r+) + ■ ■ ■ 



fyo) 

+ 9^(1) ( r ~ r +) + 



xy(0) 

(r-r+ 



(iJ) 



(5.7) 



Regularity is seen to imply that all the constants of integration are fixed in terms of 
W(o)^ a y(o)'9xy(o)}- In Particular, the g^ are not independent. 
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The general asymptotic form of the solutions as r — > oo is found to be 5 

^ - r + r 2 + ffxj/ - ^(O) + + r 3 + 

A (iJ) 

For the bulk pseudoscalar </?, the absence of a source requires that we impose <f> ^ = for a 
A = 2 operator or = ^ f° r a A = 1 operator. We wish to source only the homogeneous 
mode of the A y component of the Maxwell field and furthermore at leading order in A only, 
giving the boundary condition A^^ = 5 i,0 S^ ,Q . 

The situation is slightly more involved with the constants Gxy(o) an< ^ ^ty(o)' ^ rs ^ 
thought would be that one needs to set all these constants equal to zero since they represent 
non-normalisable modes for the bulk fields. A simple counting, however, reveals that this 
would overdefine the boundary value problem. At the horizon we saw that fl^*^ was fixed 
by regularity in terms of the other integration constants. Consequently, at the boundary, 
the constants G^'^ will be uniquely fixed after specifying G^j^ , and demanding regularity 
at the horizon, and will in general be nonzero. 

We will instead choose the constants G^ y ^ and G^ y ^ such that after a large coordinate 
transformation we will only have a source for the current. More specifically we impose 

jkG$ 0) -tu,<%$ 0) = 0. (5.9) 

We then perform the coordinate transformation 

y -> y + e~ tujt ^ Xh^ cos (jkx) , (5.10) 

hi 

where h^'^ are functions of r only. These functions are required to satisfy: 
• Near the AdS^ boundary they have the fall off 

h (hi) ~ J-G (ij l + O (r- 5 lnr) . (5.11) 

j k x y (°) v ' 



Near the horizon r = r + they have the infalling behavior 

l (0) + n (i) 



h (i,j) = f-£r ( r - r+ ) hffi + /iji'f (r - r+) + ■ • ■ . (5.12) 



The fall off condition in the UV: 

5 We chose to have m 2 — —4 in the action (4.1). Reinstating the asymptotic AdSn radius L, this 
corresponds to a mass squared m 2 L 2 = — 2. 
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• Removes the non-normalisable terms from the gt y and g xy components of the metric 
provided we satisfy (5.9). 

• It doesn't change the values of the VEVs for the various fields appearing in our ansatz. 

• It doesn't introduce new deformations coming from the g r i components. 

A particularly natural example of a coordinate transformation satisfying all of the above 
properties is one that maps our perturbations into the de Donder-Lorentz gauge 

V M <5^ = 0, V^ = 0. (5.13) 

Here V is the connection with respect to the background solution we presented in §5.1. The 
gauge conditions (5.13) become second order differential equations for the functions h^'^. 
The resulting solution satisfies all of our above requirements. 

The quantity we are interested in calculating is the optical conductivity (3.1). After 
taking into account the source of the gauge field specified by AyjjQ = 5 % '°S^ , the retarded 
Green's function function is given by (up to a constant of proportionality) 

G% Jy (^k = o) = ^^4S- ( 5 - 14 ) 

As we will explain in more detail in §5.3, for our purposes it will be sufficient to just consider 
the second order correction in A 

Gjyjy («, k = 0) = 4°'°> + A 2 A™ . (5.15) 

The first term A y ^ is the, by now well understood, current-current correlator evaluated in 
the extremal electric Reissner-Nordstrom background in Einstein-Maxwell theory [40] . 

Our aim is therefore to solve for the function a^ 2,0 ' 1 and more specifically we would like to 
extract the constant from it. After staring at the equations of motion one will notice 

(2 2) (2 2) (2 2) 

that in order to achieve that we will not need to solve for the functions a y , g\ y ' , gxy and 
(p( 2 > 2 ) . The eight functions we will be solving for are ay°'°\ afy l,1 \ o, y 2 '°\ ip^ 1,1 ^, 9ty'°\ 9ty ^ '> 
gty'°^ and g&jj 1 ^ which satisfy five second order differential equations and three first order 
ones. In order to find a solution we have five constants of integration from the horizon data 
{ a sKo)^ ' a w(oV ' "s^o)^ ' ^(o) 1 ^ ' ^y'^o) } and eigllt constants of integration from the asymptotic 
boundary {a™ , A™, A™,*$\ ^'g), Ggg> , G^g, egg} }. This counting reveals 
that we will now have a unique solution for the perturbation once the background and the 
source have been fixed, including the frequency uj. 
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5.3 Numerical results 

Having set up the problem we would like to solve, we now employ a numerical shooting 
method appropriate for the case at hand. As we have outlined in §5.1, the first step of 
the calculation is to correct the background up to second order in A. To achieve this, we 
solve the equations of motion perturbatively order by order in A which, as we explained, 
boils down to solving a system of six second order and one first order differential equations. 
The background solution will be uniquely specified by a set of thirteen constants c\ = 

\^u(0)' S ' ^aa(O) ' °(0) ' a (0) J anCL ° 2 ~ \ U xx{Z) ' U n(3) ' ^(1) '^(1) J" 

In order to specify these constants, we numerically integrate the equations of motion 
from a distance close to the horizon r = r+ + e up to a finite radius in the bulk r = r m , 
using the near horizon expansion (5.4) as the initial conditions at r = r + + e. Following this 
procedure allows us to express the functions and their derivatives at r = r m as a function 
of the constants of integration ci . We also integrate numerically the equations from a large 
radius r = R down to the same point r = r m using the asymptotic expansion (5.5), this time 
expressing the functions and their derivatives as a function of the constants of integration 
C2- In order to get a solution in the full bulk we need to match the values of the functions 
and their derivatives we have found at r = r m by suitably adjusting the constants c\ and 
C2- A simple counting shows that this solution for the constants c\ and C2 will be unique 
i.e. the number of equations we will have from matching the two solutions will be thirteen. 

Once we have specified the background we follow the same procedure to solve for the 
perturbation we described in §5.2. This time we will be using the near horizon expansion 
(5.7) and the asymptotic expansion (5.8). After matching the two solutions in the middle of 
the bulk, as we did for the background solution, we will have uniquely fixed the other thirteen 
constants of integration: the background constants ja^y*, a tfp) ' ^(0)^' ^it/'(o)} anc ^ 

the fluctuation constants A§ A™, *$\ G^, 0^,0^,0^. In total 

we will have to specify twenty six constants of integration, thirteen for the background and 
thirteen for the perturbation. The constant we are ultimately interested in is A™ which 
will be precisely the A 2 correction of the current-current two-point function, according to 
(5.15). The A 2 correction to the conductivity will then be given by 

^(2,0) 

CT (2) = ^«_ (516) 

From the spectral analysis of §4.1 and the matching argument, we expect that if an IR 
relevant or unstable operator mixes in the perturbation, we should see its imprint as a strong 
contribution to the d.c. conductivity as we lower the system's temperature approaching 
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extremality. To illustrate this we have fixed n = 36 and c\ ~ 8.47 in the Lagrangian 
density (4.1), recall we already put m 2 s = —4. This choice gives exactly one value of 
momentum k c ~ 1.27 for which the corresponding IR mode saturates the AdS2 BF bound, 
and no unstable modes. 6 The explicit expressions for the spectrum around the AdS2 x R 2 
background solution is more complicated than the case with n = 4 that we chose for 
illustrative purposes in section §4.1. In figure 4 we present the dominant v{k) for this 
specific case demonstrating that at precisely k c we have a marginally stable mode with 
v = 0. This tuning to the boundary of instability is one way to address the fact, discussed 

0.6 
0.5 
0.4 

v 

0.3 
0.2 
0.1 
0.0 

0.0 0.5 1.0 1.5 2.0 

Figure 4: Plot of the dominant exponent v as a function of k for n = 36, m? s = —4 and 
c\ ~ 8.47 illustrating the existence of a mode with v = at k c ~ 1.27, and no unstable 
modes. 

at the end of §4.1, that the model we are studying typically has an instability, with the 

effect of dialing ki/Vj* being to determine whether the instability manifests itself in the 

conductivity or not. Though we will quote values of ki in this section, as measured in the 

IR theory which as z = oo and hence dimensionless momenta, we should recall from §4.1 

that the invariant UV quantity is ki/\^P, that can be tuned via doping. Ultimately one 

would like to extend the model so that the system can be tuned through the instability. 

Fixing a small value for uj (we take Ld/fJ, = 10 -6 ) and solving the equations of motion 

as we described above, we have calculated the correction to the d.c. conductivity due to 

lattice induced scatterings and we have plotted the result against temperature in figure 5. 

As we have explained in previous sections, the expectation is that at low temperatures this 

should be dominated by a power law T 2v . The behaviour we observe for the two 

6 We have checked that for these values of the coefficients, there is also no instability localized away from 
the near horizon region. 
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different values of lattice momentum we chose to study precisely matches our expectations. 
In particular, the critical case gives a ~ 1/T. We see that quite low temperatures are 

0.0 

-0.2 

-0.4 

dlog(Re(o-)) 

diog(r) _ 06 

-0.8 
-1.0 

0.00 0.01 0.02 0.03 0.04 

Figure 5: Plot of the logarithmic derivative of the temperature dependence of the 
d.c. conductivity for n = 36 and c\ rj 8.47. The top (red) curve has k^ = 1.5 and the 
bottom (blue) curve has ki = 1.27. These correspond to v = and v ~ 0.2, respectively. 
At low temperatures the expected scaling <j( 2 ) ~ j^f— l j s recovered. 

necessary in this model to reach the true scaling regime. These lower temperatures are 
costly to reach numerically. Figure 5 confirms that the various genericity assumptions we 
have made in the matching arguments are legitimate. 

We now turn to the correction to the optical conductivity. For the same values of the 
constants n and c\ as before, we have fixed the lattice momentum to be precisely the value 
for which we have the marginally stable IR mode participating in the perturbation. We 
have plotted the correction of the optical conductivity as a function of the frequency u in 
figure 6 for two different low values of temperature T. We present these results in two plots. 
The first shows the optical conductivity up to high frequencies. The features one notices 
are an extended tail at low values of the frequency u and a regime of negative values. The 
negative regime implies that spectral weight is being transferred at interband energy scales. 
We will see in the following paragraph however that most of the spectral weight in the broad 
tail at intermediate energy scales is being extracted from the Drude peak. The second plot 
shows that the broad tail exhibits a power law falloff, again in agreement with our general 
expectations and matching arguments. If we fit the power law, we do not quite find the 
expected a ~ 1/u for this critical hi,. This is because we are not at low enough temperatures 
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to have reached the true scaling regime, as we saw in figure 5. We have checked that as 
the temperature is lowered, the scaling tends towards the anticipated a ~ 1/oj. We can 




Figure 6: The correction to the real part of the conductivity Re (o^ 2 **) for n = 36, ml = —4 
and ci ~ 8.47 at = 1.27. The red (lower) curve has T//i ~ 0.01 and the blue (top) curve 
has T//i ~ 0.006. The left plot shows the redistribution of spectral weight at interband 
frequencies, while the right, log-log, plot shows the expected intermediate-low frequency 
scaling regime. 

recall that the zeroth order optical conductivity (without a lattice), to which the results 
of figure 6 should be added, shows a featureless increase from zero at low frequencies to 
a constant at frequencies above w [40]. In fact, it is a general feature of known extremal 
holographic backgrounds that the optical conductivity vanishes as a positive power of u at 
low frequencies [43]. Additional dynamics, such as the mixing with a critical mode we have 
implemented here, seems to be necessary to obtain an extended quantum critical tail such 
as that in figure 6. We once again emphasize that this tail is not a Drude peak, there is 
still a delta function at lv = 0, as we now discuss. 

Let us examine the spectral weight transfer more closely. The conductivity sum rule 
implies that the integral of the real part of the conductivity over the positive real axis should 
not change upon introducing our ionic lattice. Therefore, at each order in the A expansion 
we can expect the correction to the conductivity to integrate to zero. This integral includes 
the contribution from the Drude peak delta function. The correction to the delta function 
can be read off from the zero frequency limit of the real part of the current current correlator. 
We can write 

RecJ (2) = -vrReG^OHM + M . (5.17) 

id 

The left hand plot in figure 7 shows the real part of the current correlator at zero frequency 
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as a function of temperature for two values of v. From this plot we see that there is delta 
function with negative weight in the correction to the conductivity. This in turn shows that 
there is spectral weight transfer from the Drude peak to the broad tail. The spectral weight 
transfer to the broad tail increases with T in one case and decreases with T in another. The 
right hand plot in figure 7 shows the integrated spectral weight 

,ImG( 2 ) {(J) 



7T 



ReG (2) (0)+ / doj' 
2 Jo+ 



UJ 1 



(5.18) 



plotted at the critical momentum for the two different temperatures considered. The func- 
tion asymptotes to zero, confirming the sum rule and providing a nontrivial check of our 
numerics. From the integrated spectral weight we see that the weight lost from the Drude 
peak already balances out the broad tail contribution, while the negative region we noted 
before in figure 6 at higher frequencies is then balanced out by higher energy processes, and 
so seems not to be closely related to the critical power law tail. Spectral weight redistribu- 
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Figure 7: Spectral weight transfer. The plots have n = 36, ml 
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Left, the zero frequency limit of the real part of the Green's function correction Gjy Jy as a 
function of temperature. The bottom (red) curve has k^ = 1.5 and the top (blue) curve has 
ki = 1.27. The negative of this curve is the spectral weight extracted from the Drude peak. 
Right, the integrated spectral weight defined in (5.18). The red curve has T j \x ~ 0.01 and 
the blue curve has T/fi ~ 0.006. Both are at the critical lattice spacing ki = 1.27. We see 
that the sum rule is satisfied after integrating up to a few times the UV scale \i. 



tion from Drude peak to broad tail is a key experimental feature in unconventional metals, 
e.g. [3, 4, 5, 44]. 

Spectral weight redistribution has recently been discussed in the context of fermionic 
quantum BKT transitions [45, 46]. There it has been suggested that the discrete scale 
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invariance associated with the complex scaling dimensions in the unstable phase may be an 
underlying principle of the spectral weight redistribution occurring at Mott transitions. 

In the remainder of this section we shall make some comments about the case in which 
there is an actual instability over some range of momenta. A first question to address 
is whether spontaneous homogeneous currents are generated via coupling to the lattice. 
According to the analysis around the Reissner-Nordstrom black hole of [30], the extremal 
near horizon version of which we repeated in §4.1, we will have unstable modes at finite 
momenta. A perturbative analysis in A reveals that these modes will be unstable around 
the backgrounds we are considering here as well. It was shown in [30] that a particular 
mode with wavenumber k c will first become unstable as we lower the temperature. Due 
to the presence of a lattice, we will now have the mixing of all modes with wavenumber 
k c + nki where ki is the lattice wavenumber and n is integer. One of the fields that 
will participate in the unstable perturbation will be the gauge field. The question now is 
what would happen in the case where k c = mkL for m a positive integer. This is the 
critical case we are most interested in. Is a spontaneous homogeneous current generated? 
More generally, this question can be asked whenever m k^ lies within the range of unstable 
momenta. We considered this case and explicitly constructed numerically the unstable 
mode at next to leading order in A. Doing so revealed that the expectation value for the 
homogeneous current is zero to within our numerical accuracy. It seems that this result 
should be provable from general principles. 

We have also studied the behavior of the conductivity close to the temperature where 
a striped phase transition occurs, assuming that our lattice is precisely such that k c = 
mki so that the unstable mode is directly communicated to the current. We find that 
the d.c. conductivity diverges as Re (<r^) oc (t — t] ^ where is specified by the 
existence of a static mode around the A = background. The presence of a lattice will shift 
the critical temperature: T c = + A 2 + • • • [47] . The full correlator should diverge 
at T c . 

6 Discussion 

In this discussion section we will reiterate the main points of our paper and simultane- 
ously contextualize our results in terms of previous holographic and other theoretical work, 
highlighting potential experimental features. 
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6.1 Broader perspective on linear resistivity, BKT and other models 

In §2 we described quantum BKT transitions from the standpoint of the near horizon 
geometry of an extremal black brane. In particular, the scaling dimension of the operator 
O took the form A = (2 + z)/(2z) + u, with v the square root of a quantity that becomes 
negative beyond the critical point. This type of formula is typical in holographic frameworks. 
The fact that 2 A = 1 + 2/z when v = underlies the scaling of the spectral weight at the 
critical point in (2.5): lim w _>o im Gqq (oj) /oj ~ 1/T. In the large ./V limit, in which the 
holographic description is classical, this implies that the 'double trace' operator OO is 
marginal. 

A more general field theoretical description of quantum BKT transitions [19] is as a 
merger and subsequent annihilation of two fixed points of the RG flow. At the transition 
point, where the two fixed points merge, one can expect in general that there should be 
a marginal operator (that does not exponentiate - this statement at the critical point is 
in addition to the marginal operator that is presumably necessary to dial the dimensions 
of operators). It seems likely therefore that OO will be marginal in general. However, 
away from the large N limit, or some other kind of generalized mean field description, the 
dimensions of operators do not add and therefore we cannot draw any conclusions about the 
dimension of O. It follows that O will not in general have the dimension A = (2 + z)/(2z) 
at the critical point, ultimately not leading to a linear resistivity. A BKT transition in 
itself therefore does not appear to be sufficient. On the other hand, if the operator that 
condenses does admit a mean field description, then our mechanism for a linear resistivity 
will automatically go through. Theories with holographic duals are one set of concrete 
examples where this is the case, even when the formula for the operator dimension is more 
complicated than (2.3). Other examples were given in [19]. Note that "mean field" here 
does not imply that the transition is of mean field nature - it is not - rather that there 
exists a description in which the order parameter may be treated as a classical variable. 

Beyond BKT transitions, previously studied models have used a conceptually similar 
approach of coupling the electric current to degrees of freedom with a 1/T spectral weight. 
For instance in [48, 49] such degrees of freedom resided on a density of impurities tuned 
to a specific value of the impurity coupling for which the model becomes solvable. The 
"Gaussianization" of the model at their critical coupling is reminiscent of our observation 
above that a mean field description of the unstable operator at the BKT transition is 
necessary to obtain the required temperature dependence of the spectral weight. The model 
also shares the feature of our approach that the Drude contribution to the conductivity is 
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hoped to be subdominant relative to the quantum critical contribution. 

A mean field limit was also tied to the emergence of a linear in temperature resistivity 
in Kondo lattice models, e.g. [50, 51, 32, 52]. These models are not dissimilar to those men- 
tioned in the previous paragraph. A 1/T spectral density of the localized spin excitations 
is fed into the self energy of the conduction fermions, leading to marginal Fermi liquid phe- 
nomenology [14]. Holographic realizations of similar models are mentioned in §6.3 below, 
see [52] for an extended comparison. 

In short, the general strategy is to find a well-motivated 1/T spectral density and then 
transfer this temperature dependence to the conductivity. The holographic BKT transitions 
we have discussed provide a robust way to this, using irrelevant operators rather than 
explicit fermions to transfer the spectral weight, with the interesting behavior occurring at 
the boundary between a critical phase and an ordered phase. 

6.2 Experimental signatures of the mechanism 

Incorporating the effects of a Drude peak will presumably be necessary to connect with key 
features of the phase diagrams of unconventional materials, such as the crossover to Fermi 
liquid behavior away from the critical point. However, it seems possible to distill at least 
three generic ingredients of the scenario we have outlined that may manifest themselves 
experimentally in one way or another. 

Firstly one could look for evidence of BKT scaling in quantities such as the critical 
temperature (2.6) or the expectation value of the condensed operator. In general, on the 
side of the critical point in the phase diagram where an operator condenses, one expects 
quantities to be controlled by a dynamically generated scale with an exponential dependence 
on the distance to the critical point [19]. 

A second generic feature is likely to be an asymmetry in the optical and d.c. conductivities, 
as well as other quantities that couple to the unstable mode, between the two sides of the 
quantum critical point. On the side where an operator condenses, at sufficiently low tem- 
peratures and frequencies, one expects only logarithmic corrections to the 1/T and 1/uj 
scalings of the critical point itself. See appendix A. The other, disordered, side of the 
critical point should be characterized by varying power law dependences along the lines of 
equations (3.6) and (3.7) above. 

Thirdly, the nature of the transition, in which a scaling dimension becomes complex as 
two critical points annihilate [19], requires the existence of a critical phase in the IR, of 
which the BKT quantum critical point is a boundary. This critical phase should manifest 
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itself in scaling behaviors such as (3.7) persisting down to arbitrarily low temperatures on 
the disordered side of the critical point. This scaling behavior may or may not require 
disentangling from other physics, such as superconductivity, present in the phase diagram. 

6.3 Context: linear resistivity in holography 

There exists a large body of interesting previous holographic work attempting to understand 
the origin of a linear in temperature d.c. resistivity. These works can broadly be classified 
into two classes, depending on the fate of the translation-invariance delta function in the 
conductivity. When the conductivity is due to probe branes [53, 54] or probe fermions 
[55, 56], the degrees of freedom of interest are taken to be a parametrically small part of a 
bigger system. In this case the charge carriers can dump their momentum elsewhere and 
it does not return on the 'experimental' timescale. A second option, that avoids diluting 
the carriers, is to explicitly break translation invariance via parametrically heavy degrees 
of freedom such as impurities [36] or a lattice [24, 34]. 

For both classes of models, the d.c. conductivity typically depends on parameters of 
the model. The dynamical critical exponent z in [54] and the dimension of IR operators 
in [55, 56, 36, 24]. These models do not provide a compelling motivation, at least without 
additional input, for selecting the exponent that results in a linear resistivity. An interesting 
example of how additional input can help is [57], where the dimension of the impurity 
operator of interest is protected by supersymmetry. It is furthermore unclear if the probe 
limit is a reasonable approximation for the real world systems of interest. 

Our model has not addressed the fate of the delta function, but rather assumed that 
the corresponding Drude peak is swamped by a quantum critical contribution to the d.c. 
conductivity. We have suggested that this may be a reasonable starting point to discuss 
some bad metallic phases of matter. In contrast to the other models we find a linear in 
temperature resistivity occurring robustly and precisely where we would like it to: at a 
quantum critical point mediating the onset of an instability. The result is universal in the 
sense that it does not depend on the specific microscopic realization of the BKT transition. 
It does require, however, the existence of an irrelevant operator that can communicate 
the instability to the current correlators. We achieved such a coupling in §4 by having a 
vectorial instability at a finite wavenumber that was was communicated to the current via 
an irrelevant lattice. 
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6.4 Open questions and future directions 

The essential feature of the holographic dynamics we have discussed is the BKT quantum 
phase transition. The realization of such a BKT quantum phase transition in e.g. some 
Hubbard-like lattice model as a function of doping would of course be very exciting in 
bridging the gap between extremal horizons and the chemistry of unconventional materials. 

Superconducting phases have been absent from our discussion thus far, and yet are a 
central feature of the most interesting bad metals. It seems reasonable to hope for a unified 
theory of bad metals and unconventional superconductivity. We can imagine two ways 
in which our mechanism could interplay with a superconducting instability. The first, a 
variant on a fairly conventional viewpoint, is simply that the superconducting instability 
could be triggered by coupling to the operator that develops the 1/T spectral weight at the 
critical point. A second possibility is that the critical phase itself, out of which the BKT 
transition emerges, is tightly entwined with the existence of a superconducting instability. 
It is common in holographic frameworks, see [58] for a recent discussion, for the IR of 
the superconducting phase to exhibit an emergent scaling symmetry. This could itself be 
the critical phase in which the BKT transition occurs. Such a picture would be, loosely 
speaking, in line with claims that the linear in T conductivity is associated to the existence 
of superconducting phases [59]. These critical issues deserve further thought. 

Clearly a pressing issue is to establish (experimentally?) whether the picture we have 
outlined in figure 1, namely that in bad metals the Drude peak is swamped by an extended 
'critical' contribution, is reasonable in some circumstances. An alternative possibility is that 
the extended tail in fact falls rapidly to zero at to ~ 0, so that for the strict d.c. conductivity 
the Drude contribution always dominates. This would incorporate naturally the fact that 
the d.c. conductivity is insensitive to the melting of an explicit Drude peak in the data. 
Against this scenario we can note that, at least in the holographic models we have studied, 
the quantum critical contribution to the conductivity does extend all the way down to 
u = 0. An alternative possibility would be to build a model in which momentum relaxational 
process were built into the critical sector that the undergoes a BKT transition. In such a 
model we might hope to have our cake and eat it, simultaneously having a linear resistivity 
and incorporating the Drude peak physics. 

Drude physics could explicitly be incorporated into our model by making the lattice 
potential (4.19) periodic in both directions. One could imagine solving such a model nu- 
merically along the lines of [34] . Because the lattice is irrelevant in the IR, the momentum 
relaxation rate is very slow and hence one's first thought is that there will always necessarily 
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be a well-defined Drude peak [24]. However, sufficiently close to the BKT critical point, 
there is a further long timescale in the problem, the lifetime of the mode that is going 
unstable, and this may be able to transfer spectral weight out of the Drude peak and into 
the critical tail. We made some preliminary remarks about this process in §5.3. A more 
comprehensive holographic study of this point is desirable and has the potential to shed 
light on the nature of bad metals. 

The particular model of a BKT transition that we constructed in §4 relied on local (z = 
oo) criticality to obtain critical physics at nonzero wavenumber that could then be coupled 
by a lattice onto the homogeneous electrical conductivity. It would be very interesting 
to see if this requirement can be relaxed by combining for instance BKT transitions with 
Fermi surface physics, or by finding a way to communicate BKT transitions to the electrical 
current without going via finite wavevector modes. 

The essential role of holography in our work has been to provide a setting where BKT 
transitions occur naturally and where the interplay between instability in a strongly cor- 
related medium, lattices and conduction can be studied in a theoretically controlled way. 
Holography is less able to incorporate more conventional Fermi liquid physics in a useful 
way. Any strategy to describe the full phase diagram and the crossover the Fermi liquid 
behavior will presumably have to involve a mix of distilled insights from holography and 
conventional Fermi liquid computations. 

Finally, the focus throughout this paper on the linear in temperature d.c. resistivity, 
and possible application to bad metals, should not obscure the more general points we have 
made. Two further take home messages are firstly that strong optical and d.c. conductivities 
can be achieved by coupling critical modes on the verge of an instability to the electrical 
current via irrelevant operators. Secondly, that (semi) local quantum criticality enables the 
effects of instabilities at nonzero wavevector to be efficiently coupled to e.g. the conductivity 
via lattice scattering. 
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A Matching for low frequency and temperature correlators 

Consider a spacetime with an IR (finite temperature) scaling geometry and that is deformed 
by irrelevant operators to flow to an AdS± UV fixed point. Perturb the spacetime by a large 
number N of fields & . These fields are taken such that they are decoupled in the IR scaling 
geometry but generically all couple in the full spacetime. That is, they are coupled by the 
irrelevant operators participating in the RG flow described by the full background. 

Solve the equations in the IR geometry, impose infalling boundary conditions, and then 
expand the solution near the boundary of the IR region (call it p — > 0). The fact that the 
fields are decoupled in the IR implies that 

&(p) ~ c'(> +2 /*- A ' + p A ' g?(u>,T)) . (A.l) 

Here Qf{uo,T) is the IR Green's function, the c 1 are constants, and Aj is the IR dimension 
of the operator. This is equation (3.2) in the main text. 

Away from the IR region, to leading order at low temperatures and frequencies we may 
set co = T = 0. We assume that the fields satisfy second order equations in the radial 
direction. These may be, for instance, gauge invariant combinations of the bulk fields. 
Therefore there are 2N independent solutions. From their behavior near the UV boundary 
these break into N non-normalizable and N normalizable solutions. We wish to consider 
the response to a single source. Therefore we consider solutions of the form 

N-l 

= A<S>i n + B$l + £ a m < . (A.2) 

m=l 

Here n is the particular non-normalizable solution we wish to source, is the conjugate 
normalizable solution and the <3?^ are the remaining normalizable solutions. A,B,a m are 
numerical coefficients. The UV Green's function we wish to compute is 

G R oc | . (A.3) 

The UV Green's function is now determined by a matching procedure given by a gen- 
eralization of e.g. [26]. This is possible at low frequencies and temperatures. Expand the 
'far' solutions of the previous paragraph in the regime where they overlap with the near 
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horizon solutions. One obtains 

~ « 7 _p 1+2/ ^ + 4p A/ > (A.4) 

^ ~ PL P 1+2/z - Al +Pi P A < , (A.5) 
<*4 ~ T^V+^+T^. (A.6) 

Here the a± , (5± , 7^_|_ are numerical coefficients that can be determined in principle by 
solving the differential equations in the far region. They do not have any temperature or 
frequency dependence. Matching with the IR data (A.l) now implies 

c 1 = Aa I + + B^ I + + Y,a m l I rn+, (A.7) 

m 

c I g?(cu,T) = Aa I _ + Bf3 I _ + J2 a rnl m -- (A.8) 

m 

There is no sum over / in this equation. Taking an appropriate difference between these 
two equations gives 

Av 1 + BW 1 + Y,a m z mI = 0, (A.9) 

m 

where 

v 1 = a I + Q?{u,T)-a I _, w 1 = (u, T) - , z mI = lm+ G?{^ T) - lm _ . (A.10) 

The matrix z mI is an (N — 1) x N matrix. It necessarily annihilates a specific vector 
p 1 . The Green's function we are looking for can therefore be written 

G^cx^cx-^. (A.11) 
A w -p 

In practice, we can find p by fixing one of its components and then determining the remain- 
der by inverting a square submatrix of z. We see that the Green's function is a ratio of 
polynomials in Qf"{u,T). 

We now wish to take the imaginary part of the Green's function. For real A/, the 
coefficients , j3± , Jm± are rea b whereas for complex A/, that is for uj pure imaginary, 
then the two coefficients appearing in each line of (A.6) are complex conjugates of each 
other. The case of real vj is simplest. In this case Qf ~ max(w, T) 2ui , which is small for 
temperatures and frequencies that are small compared to the UV scale. At a sufficiently 
generic point, in particular, away from instabilities that are not localized in the IR geometry, 
none of the a±,/3±, r ) m ± vanish. Then one can simply expand the ratio of polynomials (A. 11) 
in the IR Green's functions. The result to leading nontrivial orders will have the form 



G R (u,T) = <P + Y J d I G?(oJ,T). (A.12) 
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Because all the constants involved in this expression are real, taking the imaginary part we 
obtain 

ImG R {uj,T) = Y,d I ImG?{oJ,T) , (A.13) 

as quoted in the equation (3.3) in the main text. In deriving this expression we could have 
allowed for the analytic oj and T dependence of the coefficients. 

When some of the vj are imaginary, the Green's functions are not small and so we 
cannot expand the ratio of polynomials in (A. 11). In the d.c. limit wCTwe can expand 
the IR Green's function in w. Because, in the stable phase T > T c , the full Green's function 
must vanish as oj — > 0, the result will take the form 

ImG%;,T) = ^(logT), (A.14) 

for some in general complicated function F. 
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